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A model Hamiltonian is introduced which considers the main features of the experimental spin 
filter situation as s-f interaction, planar geometry and the strong external electric field. The proposed 
many-body model can be solved analytically and exactly using Green functions. 

The spin polarization of the field-emitted electrons is expressed in terms of spin-flip probabilities, 
which on their part are put down to the exactly known dynamic quantities of the system. 

The calculated electron spin polarization shows remarkable dependencies on the electron velocity 
perpendicular to the emitting plane and the strength of s-f coupling. Experimentally observed 
polarization values of about 90% are well understood within the framework of the proposed model. 
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I. INTRODUCTION 

Because of their extraordinary magnetic, optic and transport properties magnetic semiconductors^ |Among them the 
europium chalcogenides, have been subject to numerous experimental and theoretical investigationsaa. Many of them 
focused on the spin filter effect of the europium chalcogenides (SFE). 

The spin filter experiment can be arranged as a field emission experiment. Here a cooled tungsten emitter is exposed 
to a strong stationary electrical field bending down the vacuum level of the electric potential outside the emitter. The 
so formed potential barrier can quantum-mechanically be penetrated by the conduction band electrons of the emitter 
metal. The probability for this tunneling process is known to depend exponentially on the barriers height. This strong 
dependence can be used to obtain a " spin filter" by covering the original emitter with a layer of a ferromagnetically 
ordered material, e.g. a ferromagnetic semiconductom as EuS that makes the barrier spin dependent. 

From an experimental point of view the generation of highly polarized electron beams became more and more 
interesting with the growing importance of the spin polarized electron spectroscopies which are nowadays a mighty 
tool in the field of analyzing magnetic properties of surfaces and-thui filmsQ^Q. The spin filter effect of the europium 
chalcogenides (SFE) allows for polarization values of about 90%Qfl^Ll. 

Examining the SFE theoretically is worth while, too. Europium chalcogenides (EuX, with X=0, S, Se, Te ) show 
highly interesting correlation effects due to the complex interplay of itinerant conduction band electrons and localized 
4-f electrons, the latter carrying a strong magnetic moment. The so called s-f model 

H — Hs + Hsf 

= J2 C+ Cjo- - ^ ^ cr,; • Si (1) 



turns out to describe the interaction of both electron groups quite successfuUyHEj. c^'' are the usual annihilation- 
^ I (creation-) operators for conduction band electrons with spin a at site Ri. The spin of a conduction band electron at 
,. site Ri is denoted by cr^ and Si represents the spin of the half filled 4f-shell at this site. (Notation is conventional.) 
' The simplest possible approximation to solve the corresponding many-body problem is the mean field approximation 

_ _ I 

ijCT i,(T 

^Y.^T,,-'-fj{S^-)5,,)4^c,, (2) 
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(z^ = +1, Z| = —1) that has often been used to discuss the spin filter experiment, tooETElllS. The mean-field decoupled 
s-f interaction term spin dependently renormalizes the one-particle energies of the free conduction band electrons Fig. 
|l|, giving them an explicit spin- and temperature dependence. 

Below Tc the conduction band splits due to the interaction of conduction band and magnetically ordering 4-f spin 
lattice into two completely spin-polarized sub-bands. This splitting is temperature dependent and at T = of the 
order of ghS « leV. 

For the spin filter experiment this would mean, that below Tc only a ='f-electrons could be emitted, as the sub-band 
for a =J,-electrons lies much higher then the a ^f-sub-band and J,-electrons "see" accordingly a much higher tunnel 
barrier. 

The mean field picture thus predicts a degree of polarization for the emitted electrons that should be very close to 
100 % for all temperatures below and 0% above. However, this does not agree very well with the experimental 
dataB. Main failure of the mean field approximation i^^rpbably the complete suppression of spin-flip processes. 

The knoW|le|dgeiabout the s-f model has multipliednliijiiZi since the first attempts of applying it to the SFE in the 
early 70-tiesJi3l3't3. Aim of the present paper is a new interpretation of the SFE based on meanwhile made progress 
that permits us to treat the many-body problem of the spin filter effect- even including the external electric field - 
exactly. 



II. ELECTRON SPIN POLARIZATION 



A. Polarization and Probabilities 



In this section it will be shown, how the central quantity of the spin filter experiment - the polarization - can be 
connected with the dynamic quantities of the system, the Green functions, which will be determined in the following 
section from the many-body model. 

The vector Spin-Polarization P of an ensemble of electrons is defined as expectation value of the Pauli spin- 
matrices. More easily to handle is its projection on the preferential direction of spins, the scalar quantity P {degree 
of polarization): 

where N-^ is the number of electrons with spin a =ti the number of electrons with spin a ^l. We now consider 
the probabilities for electrons with original spin a to flip their spin, i.e. the flip-probabilities pa- If for example the 
originally prepared state was |, then p^ gives the probability to measure an electron in the t-state. 
Assuming an unpolarized beam of electrons (Nq^ = -^oi) passing through a spin filter box we obtain very easily the 
polarization of the out-coming beam using (H) and the flip-probabilities Pa which are completely determined by the 
physical properties of the spin filter box: 

Nn 

= iVoT(l -Pt) + ^oiPi = — (1 +Pi -Pt): (4) 

Nn 

Ni = iVoi(l - Pi) + Nom = - (Pl -Pt)) (5) 

^(1+Pe//-(1-Pe//)) 

P = ^7 = P<^ff (6) 



dcf 

Here we introduced Peff = Pi ~Pi' the effective fiip-ratio of the spin filter. 

Asking for the polarization, we should thus try to get the flip-probabilities from a theoretical model, i.e. Green 
fimctions. 



B. Probabilities and Green Functions 



Let us consider the following example: 
What is the (non flip) probability p{t) to measure at time t an electron with wave- vector k, if at t = an electron in 
this state had been prepared? 
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The answer is givenlij by the overlap of initial state cj^(0)|0) and final state c^{t)\0). {t) (creates) annihilates 
an electron at time t with wave- vector k. The state |0) is the electron- and magnon- vacuum. 

p(i) = |(0|ckW c+(0)|0)|'. (7) 

One sees the similarity of the probability p{t) with a well known entity of many-body theory, the time-dependent 
spectral density of the one-electron Green function: 

Sk{t,0)^^{Mt),c+m+) (8) 

In our special case (T=0, n=0), the average denoted by the ()-brackets has to be taken with the |0)-state, i.e. with 
the electron- and magnon- vacuum. Therefor one of the terms in the anticommutator does not contribute. 

Spectral density (H) and Green function G\^a{E) = ((cko-; cj^.))^ are closely related. Both will be determined in 
section (0). 

From (^)and ^ we obtain immediately 

p-(i)=47r2 |5kW|'. (9) 



III. THEORETICAL MODEL 

Aim of this section is to develop a many-body model of the spin filter which incorporates the main features of the 
experimental situation. 

We make the following two assumptions, which allow us to treat the corresponding many-body problem exactly: 

1. The tunneling processes of the emitted electrons are independent from each other; i.e. the conduction band is 
nearly empty (n=0). 

2. The 4f-lattice is ferromagnetically saturated, i.e. low temperatures (T=0). 



A. Geometry 



Our calculations are done for a planar spin filter, i.e. a sandwich-like batch of mono layers, which are translational 
invariant parallel to their plane. This matches the experimental situation (field emission through a EuS-layer) and 
provides at the same time the proper symmetry to include the influence of the strong electric field. 

The treatment of a planar system is based on the decomposition of the whole system into n equivalent two- 
dimensional sub-lattices (atomic layers) with Ns lattice points each: 

Ria = Ri + r„. (10) 

Ri and are perpendicular to each other, points to the atomic layer with index a, whereas Ri points towards 
the lattice point with index i inside this layer. 

From the lattice vector Ri^ this subscription is carried on to all operators and derived quantities; greek indices 
generally referring to the layer, latin indices to lattice points within this layer: 

o.^o.„, o., ^ar/, o,k,~^oi^;^. (ii) 

Fourier transformation into k-space is reasonably defined only within the two-dimensional sub-lattices, which are 
still invariant under translation: 

0,„ = -^^e^'^^'Oko, (12) 
Oka = -i=^e-'''^'0,„. (13) 
k means through all following considerations a vector of the two-dimensional Brillouin-zonc (BZ) of one layer. 
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B. s-f interaction 



It is convenient to write (0) in a for our purposes more suitable form: 

Hsf — — - J ''^{ZaSi^Uiaa- + SiaC^a-cr^iac)- (14) 

Here the spin-operators for the conduction electrons were expressed in terms of the creation- and annihilation- oper- 
ators c\^l: 

^ '^ia = '^ta]-'^tah O-ja = C+ |Cja| (15) 

^'^L = ^^Z^niaa, (16) 

(T 

with Z| = +1, Z| = —1 and the identity for the ladder operators 
has been used. 

This representation already indicates some physics: Thcr^rst term in (^4|) describes the interaction of the z- 
components of the spins similar to the well known Ising-modelE3. Accordingly it will be called Ising-term. The second 
term in ( p^ is non-diagonal in spin-indices and thus responsible for spin-flip-processes; it will be referred to as spin-flip 
term. _ 

Interactions among the strongly localized 4-f spins might be taken into account via a Heisenberg-termO, but will 
here be neglected, as we are especially interested in the situation of conduction band electrons. Magnon energies are 
two to three orders of magnitude smaller than the other energy scales of the system as bandwidth or s-f coupling. 

For similar reasons the interaction among different conduction band electrons will not be part of the rnpdel: eu- 
ropium chalcogenides are magnetic semiconductors, ordering at fairly low temperatures (T^^"^ = 16.57 KEI). In the 
interesting range of temperatures, the conduction band is nearly empty, and thus the contribution of electron-electron- 
interaction to the total energy of the system will be negligible small. 



C. External Electric Field 



The Hamiltonian of an external electric field 

i/v = -P • F (18) 
is a one-particle-operator. Rewriting it in terms of second quantization: 

n 

= E ^'^tf C+ „C,;3., (19) 

a/} 

we have to determine the matrix elements M"j^ . h!y^ acts in the one-particle-Hilbert-space of the n-th electron; | ia ) 
are elements of a complete orthonormal basis, for instance Wanmer-states. 

Calculating the M"^^ , we have to specify the vector of electric field F, and the operator of the electric dipole 

momentum P. The electric field is assumed to act homogeneously along the z-axis, i.e. perpendicular to the EuX-film 
(see Fig. |): 

F = (Fx, Fx, F,) = (0, 0, -f = const.). (20) 

With 

p = _e^?("), (21) 
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we evaluate the matrix elements: 

M:f = -ef 0-/3|(r).|ia) 
= -e/(Ria)z 

= -efaaoStf ^^auSff. (22) 

Here we assumed the Wannier functions to be eigen-functions of the position-space operator, what should hold in good 
approximation. Additionally we used (Ria)z = ItqI — aao: the z-component of lattice vector Rj^ equals the product 
of layer index a and lattice constant oq. Furthermore we introduced the interlayer potential difference u — e/ao- 

For the Hamiltonian of a homogeneous field, acting along the z-axis, i.e. perpendicular to spin filtering EuX-film 
we eventually obtain 

iacr 

Finally we thus propose the following model Hamiltonian for the spin filter: 

H = Hs + Hsf+Hv. (24) 



IV. SOLUTION OF THE MANY-BODY PROBLEM 



A. s-f model for T = 0, n = 



The theoretical model developed in the previous section will now be solved. We consider a single test-electron 
in the otherwise empty conduction band and a ferromagnetic saturated 4f-spin system. This special case of the s-f 
model is not only fundamental foiiliu; understanding of the spin filtpE. effect as it meets the experimental conditions; 
it moreover can be solved exactlyE3Ej, even for the planar geometryEj. 

The retarded one electron anti-commutator-Green-Function 

<.(^) = ((^»-4/3.))^- (25) 

will give us among other interesting information the time-dependent spectral density, which is needed for calculating 
the flip-probabilities (||) and thus the polarization of the emitted electrons (|^): 

1 f+°^ 

Sfjt) = ^J e"^""* ^kf (^) (26) 

S^^{E) = ---s{GfJE + ^0+)} (27) 

TT 

To determine this Green function, one has to solve its equation of motion: 

EGtfAE) = n c+ J + ) + if]-; c+ .))^ , (28) 

The higher Green functions, appearing on the rhs of ( [2^ ) 



^rZAE)-{{SLck,.;c+^^))^, (29) 
KfjE)) = ((^r/cfc,-.; c+ .))^ (30) 

are to be calculated by solving their equations of motion. This procedure usually leads to an infinite hierarchy of 
coupled Green functions and their equations of motion. To get any solution, physically reasonable decouplings are 
needed. One possible treatment would be the mean field approximation where ( |29| ) can be expressed in terms of ( [25| ) 
and (|30| ) is suppressed completely. 

After doing the two-dimensional Fourier transformation parallel to the planes, the equation of motion (^) of the 
one-electron Green function G°'f^{E) reads as follows: 



5 



6 

J 



(31) 



In the case of an empty conduction band (n=0) and ferromagnetic saturation of the f-spin lattice (T=0) which 
is considered here, we can solve this exactly. Writing down the equation of motion of the so called flip-function 
f^fc jj(^)^(^0|) one recognizes, that all higher Green functions may be expressed in terms of already known ones ore 
I. We find the following expression: 



vams. 



/Ns ^ '•''i^ 
q 



G 



Q/3 



(32) 



that expresses the flip-function completely in terms of the one-electron Green function G'^^{E). Here we introduced 
the complex propagator Ba{E): 



(33) 



( ^ ) and (^ij) form a coupled system of equations. The equation of motion (^ij) is at this point obviously equivalent 
to a matrix-multiplication which may be written in the following compact formulation: 



(E — Hkcr) • Gka 



hi. 



(34) 



Hkcr is the effective Hamiltonian of the one-dimensional problem of an atomic chain perpendicular to the translational 
invariant layers. In site representation it reads 



H 



kcr — 



ai 


bi2 


bi3 


bin 


b2i 


a2 


b23 ■ 


■ ■ b2n 


&31 


&32 


as 


bsn 


b„i 


bn2 


b„3 ■ 


■ ■ an 



(35) 



with its matrix elements: 



a„ =6„„(k) + S«(E) 

baf3 = eQ/3(k) 



containing the complex, local and spin-depending self-energy 



\ jhS 



JhBa (E) 



(36) 
(37) 



(38) 



the complex propagator Ba{E) (|^ and the Fourier-transformed hopping- integrals ea/3(k). It should be noticed, that 
the self-energy for a =t-electrons is trivial. This is due to the fact, that these electrons cannot participate in the spin 
flip processes, so that their interaction with the ferromagnetically saturated 4-f spins is restricted to a rigid energy 
shift of the a =|-conduction band. 

The propagator Bq, {E) and the self-energy are layer dependent (index a) but independent of the in-plane wave 
vector k. The latter is due to our neglect of magnon energies which is not necessary but convenient and reasonable 
since they are two to three orders of magnitudes smaller than other typical energy scales of the system (bandwidth, 
s-f coupling). However, our model might be solved exactly with full consideration of the magnon energies. 

From ( |34| ) we get finally the exact solution of the many-body-system by matrix-inversion: 



h (E - Hk,.) 



(39) 
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B. Solution with Electric Field 



The Haniiltonian of the electric field ( |23| ) is a single-particle-operator: 



Expressed in terms of second quantization, i.e. by means of the electron creation- and annihilation-operators the 
electric-field-operator Hy is therefore of the same structure as the operator of the kinetic energy of electrons Hs- 
They can be summed up: 

H = Hs + Hy + Hsf =J2 T^^4o.aC,Pa + Hsf (41) 

ij(T 

The renormalisation of the hopping matrix elements by the external field is equivalent to a displacement of the 
Bloch-band centers of gravity by a position dependent amount for each layer. 
An external electric field thus renormalizes the layer dependent self-energy: 

s;^ — >E;^=^E;^-aM . (43) 

However, the principal structure of the many-body problem remains unaffected. Its solution is given by 

Gk. - fi (E - Hk.)-^ (44) 

with Hkcr defined analogously to ( |3^ ) with the renormalizcd self-energy (^). 

The matrix inversion in (E3) is done numerically for a sc-film with layers parallel to the (lOO)-face. 



V. RESULTS 



As stated above, the self-energy for a =t-electrons is trivial. This is due to the fact, that this electron group cannot 
participate in the spin flip processes, so that its interaction with the ferromagnetically saturated 4-f spins is restricted 
to rigid energy shift of the a =t-conduction band. For that reason we now focus on the discussion of the a =|-results. 



A. Spectral Densities 



With the Green functions given by (44) we can calculate the spectral densities S'^". They depend on the spin 



direction a of the test-electron as well as on the layer-index a and the in-plane wave- vector k: 

S^^{E)^--^{Gi''^{E + zQ+)} , (45) 

TT 

where GJJ" is the element {G-^^)aa of the Green function matrix (^). 

Fig. 1^ shows the spectral density S^^^^ for the center layer (a = 3) of a five-layer-film with and without (right/left 
column) external electric field for no, medium and strong (top/ middle/ bottom line) interaction between 4f-spins 
and conduction band. The top left picture of Fig. ^ shows accordingly the well known dispersion of the free Bloch 
electron gas in the two-dimensional Brillouin zone. 

As the system consists of five layers one might expect five excitation energies. However, for symmetry reasons two 
pairs of them are degenerated, we see only three peaks at each k-point. 

By switching on an external electric field perpendicular to the film, this degeneracy is lifted and five different peaks 
can be observed (top right picture of Fig. |) . 

Middle and bottom line of Fig. ^ show the situation for finite values of the coupling strength J between 4f-spins 
and conduction band electrons. In each of these pictures we find two main structures: an energetically lower lying, 
broad and dispersionless so called flip band and at higher energies a sharply peaked structure which corresponds to 
the magnetic polaron and shows a Bloch like dispersion. 
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The flip band originates from spin flip processes of the conduction band electrons. Via the s-f interaction, the u =| 
conduction band electron can emit a magnon, i.e. it causes a spin deviation in the ferromagnetically saturated lattice 
of 4f-spins. As this process is of course forbidden for a =t electrons, the spectral densities S\^a] for a =t electrons 
are rather trivial. They are not shown separately. 

The sharply peaked structure corresponds to a quasiparticle with in most cases infinite lifetime: the magnetic 
polaron. This means, in analogy to the "normal" polaron, an electron, renormalized by a cloud of virtual excitations, 
namely magnons. As this quasiparticle propagates freely through the crystal lattice,, the corresponding structures in 
the spectral densities show therefore a Bloch like dispersion. 

In the strong coupling regime (J >0.2 eV) both structures are energetically clearly separated, the polaron is 
represented by (5-peaks and has therefore indeed infinite lifetime. This changes at weaker couplings: At J « 0.2eV, 
i.e. realistic values for EuS, the polaron peak touches the flip band and becomes considerably broadened near to the 
F-point of the two dimensional Brillouin zone. The lifetime of the polaron under these circumstances is finite and of 
the order of 1 S/eV « 10~^^ sec. 

The external electric field (right column of Fig. |4|) increases mainly the energetic distances between the layers and 
thereby infiuences strongly the shape of all spectral densities. The polaronic peaks become more separated. As the 
fiip band reproduces roughly the shape of the (u =t)-DOS, the flip band changes with the (cr =|)-DOS under the 
influence of the electric field, too. Because of (||) one should expect a similarly significant field induced change of the 
transition- and spin fiip- probabilities. This will be investigated in the following sub-section. 



B. Probabilities 

We apply (Q) on the typical a =|-spectral density of the s-f model. As calculated (^, U) and shown above Fig. 
^ it consists at each k-point of the broad flip band, written now as f{E), a continuous function of energy without 
singularities, and the (5- like polaron band hS{E — Eq): 

S{E) = ai f{E) + a2 hd{E - Eq) . (46) 

The a„ are the spectral weights {ai +0:2 — !)■ 
From (^) one flnds with (^ 

p(i) = |«i /(i)+«2 e-^^«f (47) 
because of the linearity of Fourier transformation. With Parsevals theorem 

\f{E)\^dE^ / \f{t)\^dt (48) 

-00 J —00 

we conclude from the normalization of f{E) that must vanish quicker than 1/t for t ^ oo. 

For a typical spectral density of our model we therefore find from ( p^ ) and (^: 

t 00 . .0 

Pit) > I a2 I (49) 

This interesting result shows, that the non-flip probability for a (a =J,)-electron is completely determined by the 
spectral weight of the polaron peak in the spectral density and therefore strongly dependent on position in the 
Brillouin zone, coupling strength and temperature. 

We ask now for the probability of finding an electron independent of its spin direction in a layer with index (3 at 
time t after we prepared it at time t = in layer a. Generalizing (|^) we find: 

pfit)=An'\sfit)\\ (50) 

We will illustrate ( ^ ) by investigating it analytically for a two layer system of free electrons (i.e. J = 0) with applied 
electric field. The spectral densities have a double peak structure: 

Si,{E) = h [ai 5{E - £;i(k)) + 03 5{E - -Balk))] , (51) 

where the a„ are again the spectral weights of the 5-peaks (ai + a2 = 1) and i?„ is the n-th excitation energy of the 
system. En is given as the n-th pole of the Green functions (Eil): 
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det(E-Hk) 



E=E„ 



£;i/2(k) = - 2e(k)-u±A 



A = 



4t2 



(52) 



A = i?i — £^2 is a measure for the energetic separation of the layers and depends on the interlayer hopping t set to 
O.leV and the interlayer- potential-difference u induced by the external electric field 
Now a simple calculation from (bOh using (|5l|) yields 



(53) 
(54) 



the probability of finding an electron in the second layer (p^^) or in the first layer (p^^) at time t after we prepared 
it at time t = in the first layer of our two layer film. Both probabilities are completely determined by A, the 
energetic separation of the layers and therefore not k-dependent. Increasing the external electric field u, we observe 
a growing confinement of thn_clectron in the layer it had been prepared in Fig. ^ This behavior is well known as 
Wannier- S t ark-lo calizat ionE£l EZl . 

Finally we want to determine the probability of finding an electron at time t in layer (3 with reversed spin — cr after 
it had been prepared at time t — in layer a with opposite spin a. These spin flip probabilities p'^^{t) are obviously 



given by the overlap of the Bloch states \tpkaa{0) ) = cj^^^|0 ) with states like |0kq„^(i) 



-q),/3, 



.it) |0). 



The spin flip process of the conduction band electron results always from a magnon emission, i.e. it is connected 
to a spin deviation in the lattice of the magnetic 4-f moments. Let us assume this magnon had been emitted in layer 
7 and carries away the momentum q. In order to get the total spin flip probability p'^^ (t) we have to sum over all 
possibilities of emitting such a (q, 7)-magnon: 



(55) 



The flip probabilities p^^{t) are related to the spectral density S^^^ of the spin flip Green function (|30|) 



J 13a 



Straight forward calculation similar to those in section 



IB) show: 



(56) 



The corresponding probability p'^^{t) vanishes, since we considered here the case of a ferromagnetically saturated 4-f 
lattice, which cannot be aligned any further and thus interdicts spin flip processes for (cr =|)-electrons. 

In a last step we have now to determine the general spin flip function -F^q " '^tiich is in the usual way connected to 
its spectral density: 



(57) 



After Fourier transformation we will get SZ^'^{t). 

The general spin flip function can be obtained from the hierarchy of the equation of motion. The calculation is 
comparable to the one in section ( [V ) and yields finally: 



^kqf (^) = ^"(^) GtliE) Gl\^{E) 



Q/3/ 



with 



The equations (|35[j3E 
termine the solution. 



dcf JhS 



Ns \JhBa{E) - 2 



(58) 



(59) 



and their generalized form which includes the external electric field (pq),(p4[) completely de- 
^) and (psh give us the probabilities for spin flip processes we were looking for. 
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The spin flip - and non flip- probabihties derived so far were evaluated numerically for a two-layer sc-film with 
layers parallel to the (lOO)-face. The results for different values of s-f coupling and strength of the external electric 
field are shown in the Figures 0j| and ^. 

It should be stressed, that the shown non-flip- and flip-probabilities were obtained in separate calculations, based on 
completely different many-body entities (one-electron Green function and spin flip function, respectively). However, 
our understanding of those complementary probabilities demands that they add up to 1 for all times and all vectors 
k (particle conservation). As the figures show, this is fulfilled in each case (bottom right pictures). 

Comparing the Figures |7[|ic| we can summarize: 

1. The spin-totalized probabilities (bottom lines) do not show a significant dependence on the strength of the 
s-f coupling J between conduction band electrons and localized 4-f moments. Accor ding ly they are equal to 
those with no coupling J = which had been determined analytically in subsection ( |VB| ) and are completely 
independent of the intra-layer momentum k. 

2. The total spin flip - and non flip- probabilities (right columns) do not show any dependence on the external 
electric field u. Accordingly they are equal to those without field u = 0. 

3. Calculations for films with various numbers of layers show: The total spin flip - and non flip- probabilities are 
also independent from the total number of layers in the film. Consequently they are the same as those which 
are calculated for a mono-layer. 

These observations simplify the practical evaluation of spin flip probabilities for the spin filter experiment essentially 
and allow to determine the polarization of the field-emitted electrons quite generally. 



C. Polarization 

According to (^ the polarization is completely determined by the effective spin flip ratio p^f f — Pi— Pi, where the 
flip probability for (cr ^D-electronspi vanishes identically and the (time dependent) fiip probability for a =:J,-electrons 
Pl is given by ( ^6|) as pointed out in the previous subsection. 

Using the observations (2) and (3) concerning pj^^ (i) made in the previous subsection the calculation of the polar- 
ization simplifies considerably: The dynamics of the spin fiip processes turned out to depend neither on the strength 
of the external electric field nor on the number of layers of the spin filter. It is completely determined by the strength 
of the s-f interaction J. 

Practically calculating the polarization therefore simply demands to evaluate the spin fiip probabilities for a mono- 
layer-film without electric field pkj,or, even simpler, its complementary quantity: the non flip probability pkj, as 

Pki +Pki = l . (60) 

In the previous subsection we showed, that after a sufficient long period of time the non flip probability pkx is given 
by the square of the spectral weight of the polaron. This "critical" time is characterized by the width of the scattering 
peak: t ^ Ih/eV — 6.6 * 10~^^sec and thus sufficiently smaller than typical amount of time an electron spends in the 
spin filter. We can therefore apply equation ( ^ ) and obtain finally: 

Pk=m = l-|«2(k)|^ (61) 

The degree of polarization is completely given in terms of the spectral weight of the magnetic polaron peak a2(k) 
and thus strongly dependent on the s-f coupling as well as the electrons intra-layer momentum k. a2(k) is obtained 
by integration of the spectral density (^). Fig. [l| shows the degree of polarization Pk calculated according to ( |6l| ) 
for different s-f interactions. 

It is interesting to see that for unrealistic small values of the coupling (J « O.OSeV), one really gets nearly 100% 
of polarization. This is due to the finite life time of the magnetic polarons under these circumstances: Polaronic 
and scattering peak touch each other near to the origin of the two-dimensional Brillouin zone, the polaronic gets 
broadened and the a =|-electrons "decay" into {a =t)-states. 

However, for realistic values of J (J^"'^ w 0.2ey) the polarization lies for all k clearly below 100% (Fig. right). 

The spin filter experiment allows only for a k-averaged measurement of the degree of polarization P. Out of all 
electrons with given total energy those electrons with maximum energy perpendicular to the tunneling barrier will be 
transmitted most likely. That are just those electrons with minimum intra-layer momentum k, i.e. the electrons close 
to the F-point of the two-dimensional Brillouin zone. According to our calculation (Fig. |ri|, right) this would yield 
P « 80%, which accounts much better for the experimental results then former mean field results do. 
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To exploit equation (^) fully, k-resolved experiments as for instance Spin Polarized Low Energy Electron Diffraction 
(SPLEED) or Spin Polarized Electron Energy Loss Spectroscopy (SPEELS) on an europium chalcogenide surface are 
suggested. Because of ( |6l| ) one would expect a significant angle- i.e. k- dependence of the degree of polarization of 
the scattered electrons (see Fig. 0). 



VI. CONCLUDING REMARKS 



Aim of the presented paper was a new interpretation of the spin filter experiments done on Europium chalcogenides 
using recent results from many-body-theory on the s-f model. This was achieved by expressing the degree of polar- 
ization of the field emitted electrons in terms of spin flip probabilities which had been determined in the framework 
of an exactly solvable many-body model of the experimental situation. This model used the s-f model in reduced 
dimensions (film geometry), including an additional term taking into account the strong external electrical field. 

Discussing the spectral densities of the one-particle Green function, the influence of the external electrical field on 
electronic behavior was seen e.g. as Wannier-Stark-ladder. 

Several physically relevant probabilities as layer resolved probability densities, fiip and non flip probabilities had been 
derived and discussed, showing for instance the Wannier-Stark-localization due to the electrical field and the strong 
dependence of the spin flip probabilities on the intra-layer momentum and the s-f coupling between ferromagnetically 
ordered 4-f spins and conduction band electrons of the spin fllter. 

The degree of polarization of the field emitted electrons evaluated in our model turned out to be well below 100% 
for all temperatures which is in good agreement with the experiments and a considerable progress with respect to 
mean field results of former works. 
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FIG. 1. Temperature dependent quasiparticle bandstructure of the s-f model in Mean-Field- Approximation (schematically) 



FIG. 2. The temperature dependent degree of polarization for field emitted electrons in a spin filter experiment. |-jThe dotted 
line is the prediction of the often r-used mean field solution for the s-f model. Data were taken from Kisker et alJj, the point 
(T=9K, P=85 %) was taken fromS. 



FIG. 3. Decomposition of a linear potential into a step-like (Vi) and saw tooth-like part. The points mark the field displaced 
centers of gravity of the conduction band in the layers, a is the layer index. 



FIG. 4. Density plot of the spectral densities Sk^i of the middle layer of a five layer spin filter film with and without 
(right/left column) external electric field for several values of coupling strength J between 4f-spins and conduction band 



FIG. 5. The non-flip probability for electrons in the spin filter is given by the square of the time dependent spectral density 
of the one-electron Green function. Shown are typical shapes of the a — da-Spectral densities of the spin filter model: 
energy-dependent (left) and time-dependent (middle), and the corresponding time-dependent non-flip probability (right). 



FIG. 6. Time dependent probabilities ^ (t) (starting at p=l) and ^ (t) (starting at p=0) of finding an electron in the first 
or second layer at time t for different values of the external field (u). Time is given in h/eV — 6.6 * 10~^® sec. Straight lines 
are time average. 



FIG. 7. Spin flip - and non flip- probabilities for a two-layer film with weak s-f coupling (J — O.OSeV) and without electric 
field. At t — & a =J,-electron had been prepared in layer one. The first line of pictures shows the time-dependent probability of 
still being in a cr =J,-state (in first layer, second layer, somewhere in the system). Second line of pictures shows the probability 
of being flipped (in first layer, second layer, somewhere in the system). Bottom line is the sum over both spin directions and 
gives thus the spin independent probability of being in the (first layer, second layer, somewhere in the system). Time is given 
in h/eV = 6.6 * lO^^^sec; the interlayer momentum k runs from (0,0) to (27r, 27r) through the two-dimensional Brillouin-Zone. 



FIG. 8. Same as Fig. ^, but with external electric field {u = 0.4 eV). The total spin flip - and non flip- probabilities (right 
column) are -compared to Fig. |^ unchanged, although the individual probabilities chariged rather drastically. The spin-totalized 
probabilities (bottom line) show the Wannier-Stark-Localization discussed in section (VB): the electron is confined to the layer 
it had been prepared in. 



FIG. 9. Same as Fig. ^, but with intermediate s-f coupling (J = 0.2 eV) and without electric field. The total spin flip - and 
non flip- probabilities (right column) show -compared to Fig. R- a strong dependence from the s-f coupling. The spin-totalized 
probabilities (bottom line) however are barely changed. 



FIG. 10. Non-flip- and flip-probabilities for a mono-layer spin filter with weak s-f coupling (J — 0.08 eV). Time is given in 
h/eV = 6.6 * 10~^^sec; the interlayer momentum k runs from (0,0) to (27r, 27r) through the two-dimensional Brillouin-Zone. 
The right picture verifles the particle conservation. 



FIG. 11. Degree of polarization Pk according to ( |6l| ) for weak (J — 0.08 eV, left) and intermediate (J = 0.2 eV, right) s-f 
interaction. The right picture shows the situation for realistic values of J in an spin filter experiment: J^"'^ « 0.2 eV following 
Wachter et al.tl. 
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